Self-organizing domain structure in a driven lattice gas 
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Instead of the homogeneous ordered particle distributions characteristic to equiUbrium systems a self- 
organizing polydomain structure is found to be stable at low temperatures in a square lattice-gas model with 
repulsive nearest neighbor interaction when the particle jumps are biased by a uniform electric field. We 
have performed Monte Carlo simulations to study the domain structure varying the temperatures, fields and 
system sizes. Revision of the traditional picture on this system raises interesting problems as discussed. 
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The introduction of driven lattice gases [Q was moti- 
vated by the demands to study non-equilibrium (open) 
systems which exhibit phase transitions and seem to be 
tractable with the improvement of dynamical methods 
using the concepts of equilibrium statistical physics. In 
these lattice-gas models the thermal jumps of interacting 
particles are biased by a uniform electric field E result- 
ing in a stationary particle transport through the sys- 
tem if periodic boundary conditions are used. These sys- 
tems exhibit many interesting features (for a review see 
1^,^). Obviously, the equilibrium behavior should be re- 
produced in the limit ^ 0. From practical point of 
view these models are able to describe the effect of driv- 
ing field on the ordering processes in superionic conduc- 
tors i|. 

Now our consideration will be restricted to a half-filled 
driven lattice gas on a square lattice with repulsive near- 
est neighbor interaction. Accepting the previous nota- 
tions the nearest neighbor interaction is considered as 
an energy unit in which both the temperature {ks = 1) 
and strength of driving field are measured . The sys- 
tem has a twofold degenerate ground state in which the 
particles form chessboard-like ordered structure. In equi- 
librium one of this structure will appear as a result of a 
sublattice ordering when cooling the system through the 
Neel temperature Tjv. The feature of the corresponding 
critical transition is well described in the literature 
The early investigations apparently confirmed the naive 
expectation, namely, the universal behavior of the order- 
ing process remains unchanged in the presence of a weak 
driving field. More precisely, dynamical mean-field anal- 
ysis 0, field theoretical investigation and Monte Carlo 
(MC) simulations Q have suggested that the Neel tem- 
perature decreases with E and the continuous transition 
becomes first order above a threshold value. These anal- 
yses have shown that the order parameter as a function 
of temperature depends on the jump rate (Metropolis or 
Kawasaki) if > 0. In both cases the transition temper- 
ature vanishes for E = 2. 

In a previous paper |^ the MC simulations were rein- 
vestigated using as large as 300 x 1500 system size for 
a fixed driving field (E — 0.4) where continuous transi- 
tion was predicted. These simulations have demonstrated 



clearly that a self-organizing polydomain structure char- 
acterizes the stationary state at low temperatures. In 
other words, the applied field prevents both the sponta- 
neous symmetry breaking and the related critical phe- 
nomena. 

In this Letter we look over the recent results and give 
further evidences to support the former picture. Our 
main goal is to illuminate the variety of phenomena and 
emphasize the joints with other areas. First we briefly 
review the results of a phenomenological model which 
clarifies the appearance of some typical lengths leading 
to a size effect. It will be shown that this size effect 
resolves the discrepancy between the MC data obtained 
for small and large system sizes. We could distinguish 
three different types of domain patterns when varying 
the temperature, field strength and system size. The do- 
main size distribution shows power law behavior in the 
self-organizing polydomain structure considered as "ther- 
modynamic limit" . Finally we summarize our main con- 
clusions. 

The visualization of particle distribution during the 
MC simulations at low temperatures makes clear that 
most of the jumps take place at the interfaces separating 
the ordered regions. Consequently, the particle transport 
is also localized along the interfaces in the presence of a 
driving field. For curved interfaces the induced current 
results in an accumulation of extra particles (holes) at 
the parts where the curvature is negative (positive) . The 
interface is driven by the field with a velocity propopor- 
tional to the charge density and perpendicular field com- 
ponent. In the mentioned phenomenological model the 
interface shape, the particle density and current along 
the interface are described by single valued functions of 
time and x coordinate Beside the mentioned phe- 
nomena the deterministic equations of motion take the 
effect of surface tension into consideration as suggested 
for equilibrium systems ||l0[| . In this model the time evo- 
lutions of the interface and charge density are separable 
if £^ = 0. Namely, the charge fluctuations along the inter- 
face are dumped by a diffusion process while the interface 
motion is controlled by the surface tension. 

Although in the presence of a field the situation is more 
complicated because of the coupling, the equations of 
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motion have trivial solutions corresponding to a neutral, 
tilted, standing, planar interfaces. Excepting the inter- 
faces parallel to the field, these solutions are found to 
be unstable against periodic perturbations in agreement 
with MC simulations . According to a linear stability 
analysis the amplitude increases exponentially with time 
if the wave-length exceeds a treshold value. The am- 
plification rate has a maximum at a given wave-length 
(A"^) proportional to \/E. Following the initial (expo- 
nential) increase of the periodic components there ap- 
pears a finger formation controlled by non-linear effects. 
In this case the typical finger width can be estimated as 
A*/2. Consequently, the growing domains are splitted 
into strips whose length depends on the thermal fluctu- 
ations. Finally the system evolves into an anisotropic, 
self-organizing domain structure at sufficiently low tem- 
peratures. In this stationary state some parts of the mov- 
ing interfaces can meet and annihilate each other while 
the accumulated charges are neutralized (in part). This 
process unites two distinct domains and generates extra 
defects into the bulk phase. All these phenomena can be 
observed clearly when displaying the particle distribution 
during the MC simulations. 

The above interfacial instability is analogous to the 
MuUins-Sekerka instability observed in crystallization 
1^ , ^ . Furthermore, a similar phenomenon is found in 
the driven lattice gas with attractive interaction for the 
interfaces separating the high- and low-density phases 
p^-p^. It is emphasized that the interfaces parallel to 
the field remain stable in the driven diffusive systems for 
both the attractive and the repulsive interactions. 

Beside the preservation of this self-organizing sta- 
tionary state the enhanced interfacial particle transport 
breaks up the monodomain structure into a polydomain 
one via a nucleation mechanism as observed in MC sim- 
ulations f7|. Due to thermal fluctuations "islands" of 
the opposite phase are formed in the homogeneous ini- 
tial state. These "islands" polarized by the above mech- 
anism elongate along the field. Using the mentioned phe- 
nomenological model we could determine the variation of 
the domain area. For simplicity we have studied a circu- 
lar domain and the charge distribution along its bound- 
ary is chosen to be a stationary solution for fixed radius. 
It is found that the area increases if the radius exceeds a 
critical value, which is proportional to 1/E. 

In accordance with the above picture we can distin- 
guish three different charcteristic lengths related to the 
interfacial phenomena; i.e., the "typical width and length 
of strips" and the critical nucleon size. The MC simula- 
tions show that Rc and A*/2 are comparable while the 
"typical strip length" (or longitudinal correlation length) 
may become significantly larger when decreasing the tem- 
perature for a fixed field {E = 0.4 in Ref. |^). 

The existence of these characteristic lengths leads to 
an obvious size effect. The self-organizing domain struc- 
ture is preserved by the interfacial instability when the 
system sizes are chosen to be much larger than these for- 
mer values. On the other hand, if the MC simulation 



is performed on a small system then one can observe an 
ordered monodomain structure at low temperatures and 
phase transition can be concluded. For example, the fi- 
nite size scaling of the order parameter can be performed 
for L < AQ and one can deduce a "critical temperature" 
T]s{E = 0.4) = 0.488. For larger L, however, the scaled 
data deviate significantly from the scaling function. 

In small system the ordered stucture evolves into its 
counterpart via the nucleation mechanism as demon- 
strated in Fig. |l| (for L = 10) where the sublattice occu- 
pation is plotted as a function of time measured in Monte 
Carlo steps per particles (MCS). 
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FIG. 1. Time dependence of sublattice occupation in MC 
simulations performed at _E = 0.4 and T = 0.48 for different 
L. 

From Figure |l| it is clear that the evolution of sublat- 
tice occupation shows different behavior for L — 120 if 
the system is started from an ordered structure. Evi- 
dently, the fluctuation of sublattice occupation (around 
1/2) decreases when choosing even larger systems. 

Beside the mono- and polydomain structures we can 
distinguish a third type of domain patterns when the 
system size equals approximately with the longitudinal 
correlation length which is generally larger then A*/2 (or 
Rc). In this case closed strips are formed with interfaces 
parallel to the field. These interfaces are not affected by 
the instability mentioned above. Due to the thermal fluc- 
tuations the interfaces move randomly, the neighboring 
ones can meet and annihilate each other. The "multi- 
strip" state can be considered as a stationary one be- 
cause the nucleation mechanism retrieves the loss caused 
by the annihilation. 

The common feature of the above situations is that the 
time average of the order parameter disappears. We have 
introduced a relaxation time to characterize the disap- 
pearance of this quantity. A series of MC simulations was 
performed varying the system size for fixed field E — 0.4 
and temperature T = 0.48. The time variation of the 
order parameter has been determined by averaging over 
150 MC runs when the system is started from one of the 
ordered structure. Following the method introduced by 
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Binder and Miiller-Krumbhaar we have evaluated the 
relaxation time for different system sizes (see Fig. ^). 
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FIG. 2. Size dependence of relaxation time in a driven sys- 
tem (open squares) and for E — (closed circles) at T = 0.48. 

For small sizes the relaxation time increases with L. 
The curve has a maximum for L w 36 and tends to a 
constant value when L oo. As a comparison with the 
equilibrium system we have repeated these calculations 
for E = and L < 40. In contrary to the driven sys- 
tem the relaxation (ergodic) time increases monotonously 
with L in the absence of drive as demonstrated. Unfortu- 
nately, the systematic analysis becomes time consuming 
for lower temperatures because of the rapid increase of 
relaxation time. 

The interfacial energy gives significant contribution to 
the total energy in the driven sytem. Consequently the 
size effect modifies the specific heat drastically as shown 
in Fig. ^. The peak increases and moves to left when 
choosing L — 16, 32 and 64. For larger systems, however, 
the appearance of interfaces decreases the peak height. 
Notice that the simulations obtained for L — 256 repro- 
duce the previous data |9|. 
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FIG. 3. Specific heat vs. temperature ai E — 0.4 for difi'er- 
ent system sizes L = 16 (open squares), 32 (closed triangle), 
64 (open circle), 128 (closed circles) and 256 (closed squares). 



For stronger field the characteristic lengths become 
shorter and the polydomain structure is more tractable. 
As a consequence the MC simulations can be performed 
on smaller systems to have reliable results. The simula- 
tions have justified the absence of long range order at low 
temperatures. Furthermore, the A singularity of specific 
heat spreads away as illustrated in Fig. ^ 
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FIG. 4. Specific heat for different fields: E = (triangles), 
0.4 (open squares), 1 (closed squares), and 1.5 (bullets). 

These data suggest that the critical behavior appears 
only in the limit E' — > 0. Unfortunately, the MC simula- 
tions are not adequate to study rigorously the crossover 
from the smooth transition to the critical one when de- 
creasing E because of the extremely large system sizes 
required to avoid the mentioned size effects. At the 
same time, the results of simulations agree quantitatively 
with the predictions of a generelized mean-field analysis 
performed at the level of 6-point approximation Jl9| for 
E>1. 

The appearance of pattern formation in the present 
model is not surprising because similar phenomena are 
found in some other open systems. The best known ex- 
ample is the Rayleigh-Benard instability |2^. The cu- 
riosity of this driven lattice gas is related to the fact 
that here the pattern formation is preserved by the en- 
hanced interfacial particle transport. That is, an interfa- 
cial phenomenon is capable to prevent the ordering pro- 
cess implying drastical modifications in the macroscopic 
features. 

The polydomain structure consisting of two phases is 
topologically equivalent to a pattern like "droplets in- 
side of droplets inside droplets ..." characteristic of the 
critical domain picture of Ising model According 
to the deterministic interfacial evolution model the ap- 
pearence of characteristic lengths leads to minimum size 
for droplets. In the real system the smaller droplets 
appear as a consequence of thermal fluctuations. The 
resultant domain structure is topologically similar to a 
quenched state roughened for a given time. To demon- 
strate it we have determined the domain size distribution 
functions in both cases. Closed circles in Fig. || show the 
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average number of domains of size s in a rectangular sys- 
tem with 500 X 1500 sites for E = 0.4 and T = 0.4. Here 
the open circles are obtained by averaging over 1000 pat- 
terns developed from a random distribution after a ther- 
malization of 300 MCS in the absence of driving field 
for a fixed temperature T — 0.4. In this former case 
there also exists a characteristic domain size increasing 
monotonously with time. 
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FIG. 5. Log-log plot of the average number of domains vs. 
their size for a self-organizing polydomain structure (closed 
circle) and for a state (open circle) during the roughening 
process. 

The domain size distribution functions are generally 
sensitive to the choice of criteria defining domains and 
interfaces. In the above comparison we used the same 
criteria to demonstrate the similarity between the two 
types of states. Prehminary results suggest that this be- 
havior is not affected by the anisotropy. Further system- 
atic analysis is required to have quantitative conclusions. 

For weak fields the interfacial phenomena including 
noises affect dominantly the long range order while it al- 
lows nearly perfect order inside the domains. In this limit 
the polydomain state can be considered as a noisy surface 
evolution problem. Such an approach raises many ques- 
tions not yet investigated in the area of surface growth 
p2| . For example, this technique may give some predic- 
tions for the longitudinal domain (or correlation) length. 
Furthermore, it would be interesting to classify the inter- 
face evolution models which are able to preserve such a 
self-organizing domain structure. 

In summary, at low temperatures the present driven 
lattice-gas model exhibits a self-organizing polydomain 
structure which is topologically similar to those appear- 
ing during the roughening process. The interfacial effects 
not only preserve the self-organizing process but they are 
able to destroy the monodomain structures. The related 
characteristic lengths are responsible for the significantly 
different features found for the small and sufficiently large 
systems. The temperature- and field-dependence of the 
self-organizing state raises many questions to be studied 
in the near future. 
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